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RESEARCH MEMORANDUM 


A SIMPLIFIED MATHEMATICAL MODEL FOR CALCULATING 
AERODYNAMIC LOADING AND DOWNWASH FOR MIDWING 
WING-FUSELAGE COMBINATIONS WITH WINGS OF 
ARBITRARY PLAN FORM 


By Martin Zlotnick and Samuel W. Robinson, Jr. 
SUMMARY 


It is shown that, for the purpose of calculating aerodynamic loads 
on the fuselage, the midwing wing-fuselage combination in subsonic flow 
can be represented by a simple system of horseshoe vortices located on 
the wing with images located inside the fuselage. Using this simplified 
mathematical model, & method for calculating the lift and longitudinal 
center of pressure on the fuselage in the presence of the wing is 
presented, 


In addition it is shown how the simplified mathematical model can 
be used for calculating the downwash behind the wing and for calculating 
the spanwise lift distribution on the wing in the presence of the 
Puselage. 


INTRODUCTION 


Mutual interference between wing and fuselage has a significant 
effect on the pitching moment of the wing-fuselage combination since 
the longitudinal distribution of the aerodynamic loading on the fuse- 
lage is altered by the presence of the wing. For configurations with 
unswept wings Multhopp (ref. 1) has developed a theoretical method for 
calculating the pitching moment on wing-fuselage combinations which gives 
good agreement with experimental results. However, for calculating the 
pitching moment on swept-wing configurations, only semiempirical methods 
such as that of reference 2 are available. 


The aerodynamic loading on the fuselage in subsonic flow is consid- 
ered in references 1 and 2 to be made up of two parts. The first part, 
which is due to the fuselage angle of attack resulting from both the fuse- 
lage geometric angle of attack and the upwash angle induced by the wing, 
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has been calculated in reference 1 for straight-wing configurations and 
modified in reference 2 for the swept-wing case. The second part, which 
is often referred to as the wing lift "carried over" by the fuselage 
(which will be referred to in this paper as the induced lift), has been 
calculated analytically for the straight-wing case by Multhopp. However, 
for swept wings Multhopp's calculation cannot be applied and no other 
theoretical method has been available. In reference 2 this component 

of loading is estimated by an empirical method. 


In order to calculate the induced lift on the fuselage in combina- 
tion with a wing of arbitrary plan form, a method is suggested in the 
present paper which is based on Lennertz's (ref. 3) theoretical work. 
Lennertz's results, which are concerned only with an unswept lifting 
line passing through the axis of an infinitely long cylindrical fuse- 
lage, are, in effect, generalized so that in subsonic flow the magnitude 
and center of pressure of the induced lift on the fuselage when combined 
with a swept lifting line, or even a lifting surface, can be calculated. 
A numerical example is given in the appendix to illustrate the method. 
The effect of finite fuselage length is estimated qualitatively from 
results of an approximate calculation of the variation with slenderness 
ratio of the induced lift on an ellipsoid of revolution combined with an 
infinite vortex. The approximate calculation is shown to give results 
which agree with results of Vandrey (ref, 4). 


In calculating the induced lift on the infinitely long fuselage, 
the wing-fuselage combination is replaced by a simple system of horse- 
shoe vortices (or doublets} on the wing, with images inside the fuselage. 
This representation may also be used, with some modifications, in calcu- 
lating the lift on the wing in the presence of the fuselage and also the 
corresponding downwash. However, for calculating the downwash and the 
lift on the wing, the simplified representation of the wing-fuselage 
combination, although it is considered to be adequate, is no longer 
rigorous, and the general applicability will depend on Capon imensa 
verification. 


SYMBOLS 
A aspect ratio (wing alone) 
a maximum radius of body of revolution 
b wing span 
C wing chord 


mean chord, S/b 


O1 
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Cy local Lift coefficient 

a major axis of ellipsoid of revolution 

Lp total lirt on the fuselage 

P pressure 

q dynamic pressure 

8 wing area (wing alone) 

に | semispan of bound leg of horseshoe vortex 
Teo free-stream velocity 

Unax "maximum longitudinal velocity on surface of body 
Xep longitudinal center of pressure 

I5 vortex strength 

y spanwise loading coefficient on wing, ce; / で 
p mass density 

Subscript: 

n particular pair of horseshoe vortices 


A prime denotes that the quantity is dimensionless with respect 
to a; an asterisk denotes that the quantity is dimensionless with 
respect to  b/2; and a bar over a symbol denotes that the quantity is 
dimensionless with respect to S. 


BASIC CONSIDERATIONS 


Midwing Configurations With Axisymmetrical Fuselages 


In the major part of the analysis the fuselege is assumed to be an 
infinite cylinder. To obtain a qualitative estimate of the effect of 
the finite fuselage length, the variation with slenderness ratio of the 
lift on an ellipsoid of revolution with a wing of infinite span will be 
calculated approximately in one of the following sections. An approxi- 
mate method for making the small correction for this effect will be 
indicated. 


RE 
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Review of results obtained by Lennertz.- Lennertz (ref. 3) has cal- 
culated the lateral and longitudinal lift distribution on the fuselage 
of an idealized wing-fuselage configuration, in which the fuselage was 
represented by an infinite circular cylinder and the wing by a vortex 
having constant spanwise circulation. The vortices trailing from the 
wing tips have images inside the cylinder, and the bound vortex ia 
extended inside the cylinder to join the trailing image vortices as 
shown in figure 1. With this configuration, the boundary condition on 
the cylinder of zero velocity normal to the surface is satisfied only 
at infinity and in the plane normal to the cylinder axis which passes 
through the vortex so that it is necessary to superimpose an additional 
potential, which Lennertz calculated. For this case the lateral lift 
distribution was obtained by considering the momentum change in a verti- 
cal plane infinitely far behind the wing, and the longitudinal lift dis- 
tribution was obtained by the use of Bernoulli's equation. 


induced lift on an infinite cylinder with a finite wing having 
constant spanwise circulation distribution.- It will be shown that 
the fuselage lift Lp and its lateral distribution dlp /dy are not 
affected by the additional potential so that only the components of Lp 
and GLr /dy due to the vortex potential need be calculated. 


Neglecting all other singularities, the lift on a longitudinal 
section of the cylinder dy (see fig. 1) due to the vortex potential 
is calculated as follows: 


The pressure p of any point on the surface of the cylinder can 
be written 


P = Po + $0|(U,. + du)” + Av? 2 


where Au, Av, and Aw ere, respectively, the longitudinal, lateral, 
and vertical components of velocity induced by the vortices. The sec- 
tion lift  dLe/dy is then written 


dat, oo 
ay Te a> n)a 


where Pù and p} are the pressures on the upper and lower surfaces, 
respectively. 
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Then 


CO 
ka eeu, | Au ax (1) 
dy =% 


Only the velocities induced by the bound vortex contribute to the lift. 


It may be noted that, since the distribution of the longitudinal 
velocities induced by the bound vortex is symmetrical about its axis at 
every section dy, the longitudinal center of pressure of the Lift due 
to the vortex potential is on the axis of the bound vortex. 


A closed expression for the integral in equation (1) may be readily 
derived as follows: Consider the rectangular path indicated by the 
dashed line in the upper right sketch in figure 1 and the cross section 
downstream at infinity shown in the lower right sketch of the same fig- 
ure. The line integral of the tangential velocity component, taken 


around the complete rectangle, Pnu ds must be Ij where the path 


a2 
b/2 


where it does not link one of the horseshoe vortices (0 < M « 


« ixl < z and must be zero 


) 


links one of the horseshoe vortices ( 


a 
b/2 
Then the desired integral in equation (1), which is the longitudinal 


portion of the complete line integral [uw ds, must equal あ Au ds 
minus the line integral along the short vertical line at infinity 


i 
n Aw(e,y,z)dz, where u refers to the upper surface and 2 refers to 
u 


CO 
the lower surface. The integral Fi Au dx is 
= QO 


f. Au dx = pau ds Jf Awe ,y,z)àz (2) 


Since 
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then 
1 
n Aw dz = fa- Pi 
u 


where Pu and Ø: are the potentials on the upper and lower surfaces, 


respectively, in the plane x = œ. When 77 < |» | < lal, 
2. \ 2. e 
Bu = の Fb tandt Na? - ye ー も anー ユ Na? - y - tan" Ne ey + 

2 cid mas y 22 
2 2 b/2 

2 abis? - y? 

tan71 さき E - + I3 tan7t 5 

y + a (3 2 , a2 
b/2 2 


and 


2 
and when Ox | | < | 一 一 | the strength of the image vortex located at 
a2 
ーー- is increased by m 
b/2 
I | as | 2. \ 2 
Pu = ~P, = [P lian 12 - tan7+ 32 に - [n + tan 1 2 y^ + 
on D a2 
yrs Yrs eee 
b/2 
bf の 2 
\ 2 esla- - Y 
puer eg = L py ftant 2 -Z 
co (B - a 
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and. 


pau as = 0 


Then, from equation (2), the lateral lift distribution of the induced 
lift can be written 


2 


bV.2 
ang ap? - y 
dy 


= pU4Up|l - = tan~+ ee (-a <y <a) (3) 
(B) _ a2 
2 
and after integration over y 
2 
= 260 ( -Ê (3) 
Le ol b E 


which are the same as the expressions obtained by Lennertz for dLp/dy 
and the total fuselage lift Lf. The components of Lp and dlp /dy 
due to the additional potential must therefore be zero. 


The longitudinal lift distribution on the cylinder calculated by 
Lennertz, which included the effect of the additional potential as well 
as the vortex potential, is slightly different from that which would be 
calculated due to the effect of the vortex potential alone. However, 
both distributions are symmetrical about the axis of the bound vortex. 
The longitudinal center of pressure of the lift due to the vortex poten- 
tial must be at the bound-vortex axis because the longitudinal velocities 
induced on the surface of the cylinder by the bound vortex are the same 
in front and in back of the bound vortex as noted before. That the 
longitudinal center of pressure of the lift due to the vortex potential 
and the additional potential must be at the axis of the bound vortex as 
the Lennertz calculations show (although he does not say this explicitly) 
can be shown as follows: 


Superposition of infinite vortices canceling the semi-infinite 
trailing vortices of figure 2(a) as shown in figure 2(b) will not change 
the longitudinal lift distribution on the cylinder since no longitudinal 
velocities are induced. It is apparent that the system is the same as 
before, with the direction of the trailing legs of the horseshoe vortices 
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reversed; therefore, the longitudinal distribution on the infinite 
cylinder must be symmetrical about the axis of the bound vortex, since, 
if it were not, the longitudinal lift distributions on the cylinders in 
figures 2(a) and 2(b) would be different. 


From the analysis in this section of the wing and cylinder combina- 
tion it can be seen that, for calculating the lift on the cylinder, its 
lateral distribution, and the longitudinal center of pressure, it is 
necessary to consider only the effect of the vortex potential. This 
Lift due to the vortex potential will be referred to hereinafter as the 
"induced lift." 


Simplified representation of the wing-fuselage combination.- In 


this paper, the wing-fuselage combination will be represented by a system 
of discrete horseshoe vortices and images, so that configurations with 
Wings of arbitrary plan form may be treated. It is necessary to super- 
impose two pairs of horseshoe vortices of the type shown in figures 1 
and 3 to obtain the vortex-image system of figure 4. The vortex-image 
system of figure 4 can be used to represent wings of arbitrary plan form 
in the manner shown in figure 5 by locating the bound vortices on the 
wing quarter-chord line or by distributing them over the wing surface in 
the same manner. Since the discrete horseshoe vortex becomes a doublet 
line as the length of the bound leg approaches zero, it is also possible 
to represent the wing by & continuous distribution of doublets. 


Figure 3 shows the larger pair of vortices with span 2(Yn + 8) and 


strength +I near the smaller pair of span 2 (Yn ー g) and strength -T. 
When the bound legs of the smaller are moved to coincide with the bound 
legs of the larger, the net vortex strength along the section where they 
coincide is zero, and the remaining sections of the larger bound legs 
form the bound legs of the desired system (see fig. 4). From equa- 

tions (3) ana (4), the lateral distribution of the induced lift and the 
total induced lift for the vortex-image system shown in figure 4 are 
given by the following expressions: 


2(yn - 5)Ya? - E 2 (ya + a)Ma2 - y* 


y^ 
Ls = (=) Ee ten! —L——— ————————— - tan 
qs ? /n (Yn ー g ) - a? (Yu 4 s)? - a? 


(5) 
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pur x ge Sa (6) 
qs € a ye - s? 


cc 
where ( a is the loading coefficient on the wing at the station y,. 
n 


Since it is known that the longitudinal center of pressure of the 
induced lift (and total lift) for each pair of horseshoe vortices and 
images is still located at the axis of the bound vortex, it is possible 
by superposition to calculate the pitching moment on the fuselage in 
the presence of a wing of arbitrary plan form if the lift distribution 
on the wing is known, although the complete longitudinal lift distribu- 
tion cannot be calculated unless the component of loading due to the 
additional potential is calculated. The total induced lift and its 
lateral lift distribution can also be calculated by superposition. The 
method is discussed in one of the following sections, and an illustra- 
tive example is given in the appendix. 


Effect of Finite Fuselage Length on the Induced Lift 


To obtain an estimate of the error involved in the assumption of 
the foregoing analysis that the fuselage is infinitely long, an approxi- 
mate calculation will be made of the effect of slenderness ratio on the 
induced lift of an ellipsoid of revolution combined with an infinite 
vortex. The limiting case of the spherical fuselage is treated first 
in the manner of Vandrey (ref. 4), and then the general case is treated. 


Calculation of the lift on a sphere-vortex combination.- The poten- 
tial of the sphere combined with an infinite vortex (fig. 6(a)) is written 


る = ーー tent Z 


+ U x [1 "T TENE REM (7) 
2 x i 2 GE 


wbere the potential for ihe sphere in the free stresm is 
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(8) 


6; = Uox ll +$ 


ICE 


and the potential for the vortex is 


E (9) 


eux x 


The lift on a section of the sphere in the plane y = Constant is 


= d ra (Pu - p。) dx (10) 


where p, and p, are the pressures on the upper and lower surfaces, 
respectively, From Bernoulli's equation, 


P = Po + ぁ | + Au)? + (v + Av)? + (w + aw)?] (11) 


Pu - Pz = 2p(u Au + w Aw) 


Where Au, Av, and Aw are the velocities induced by the vortex on 
the upper half of the sphere and u, v, and w are the local velocities 
on the surface of the sphere. Thus 


2x ge . y? 2x g2 . y2 
(12) 
- 
2 € a2 su" - 
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The lift on the section at y = Constant is 


a2 


which may be written in the form 


since 
Aus 回 
2X Re ャ ー 
and since 
u = Unax 
when x= O, so that 
3 
Umax = 5 Vo 


LI 


(13) 


(iH) 


For an infinite cylinder with an infinite vortex, equation (14) will 


also hold. In this case 


Unax = YU, 
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From equation (13) or (14) the lift obtained for the sphere is 


dy n 
and for the cylinder 
ay Peb 


Equation (14) gives the exact result for the case of the sphere 
with an infinite vortex and the same result as the method of Lennertz 
for the case of an infinite cylinder with an infinite vortex. It is 
suggested, therefore, that equation (14) may be used to obtain a quali- 
tative estimate for the intermediate case of an ellipsoid of revolution 
having & slenderness ratio between l and infinity. This assumption will 
be made in the following section and it will be shown that the value 
obtained for the induced lift by using equation (14) is very close to 
that obtained by the more accurate mathematical treatment of Vandrey 
(ref. 4) for the case of an ellipsoid having a slenderness ratio of 5. 


Calculation of the lift on a combination of an ellipsoid with an 
infinite vortex.- From equation (14) the local lift on the ellipsoid is 


written 


AL, (y) 
Eom f. "LL (15) 


From figure 6(b), 


f(y) = ayi - (zy 


(16) 
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The integration indicated in equation (1h) yields 


1 ale  UmaxlU. (17) 
UT, dy | a 
PUSH V 117 


The values of Unax[ Uc given as a function of the slenderness ratio in 
reference 5 are shown in figure 7. 


Figure 8 shows a plot of — e against d/a. The result of 
P cd 
a calculation by Vandrey for the case of an ellipsoid with 2-5 is 


shown to fall near the curve in figure 8. Since Vandrey's calculations 
appear to be accurate to 0.05, the agreement may be even better than is 
indicated in figure 8. Vandrey's results are obtained by difficult 
computations for each body separately, and his method does not yield a 
general expression similar to equation (17). 


Figure 8 indicates that the induced lift on the ellipsoids having 
high slenderness ratlos (d/a > 5) is about 90 percent of the induced 
lift on the infinite cylinder. Since the fuselage is similar effec- 
tively to a semi-infinite cylinder because of the wake which extends 
behind it, the loss in lift due to the finite fuselage length is about 
half of that indicated by the calculation for the ellipsoid of revolu- 
tion. The value of induced lift obtained by assuming the fuselage to 


L Unaxl Us 
be an infinite cylinder may be multiplied by the factor 5 l + —À 
313 
to correct for the finite fuselage length; however, this approximate 
correction can usually be neglected since it is nearly unity for gt 
practical slenderness ratios. Since the correction is small, it is 
assumed that it may be applied directly for the finite-wing case without 


d 
be multiplied by equations (5) ana (6) to correct for finite fuselage 


length. 


1 Unax! Us 
introducing significant error, so that the factor を [+ + PEA is to 


1h 
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APPLICATION OF SIMPLIFIED MATHEMATICAL MODEL TO CAICULATION 


OF THE AERODYNAMIC LOADING AND DOWNWASH 


In the following sections, the method for calculating the induced 


lift will be described and discussed 
spanwise loading on the wing and the 
methods for calculating the downwash 
are not rigorous since the effect of 
has not yet been calculated, must be 
tion. 


and methods for calculating the 
downwagh will be outlined. The 
and spanwise lift distribution 
the additional potential, which 
approximated by a simple correc- 


Although the validity of the correction must depend on experi- 


mental verification, 1t ls believed to be adequate, and the exact value 
of the additional potential may be incorporated into the method immedi- 


ately when it is calculated. 


Method for Calculating the Induced Lift on the Fuselage 


If the lift distribution on the wing in the presence of the fuse- 
lage is known, say from reference 6 or the method outlined in the 
following section (see fig. 9), the induced lift may be calculated very 


simply. 


A sample numerical calculation is shown in the appendix. 


It is 


only necessary to substitute into equations (5) ana (6) the value of 


loading coefficient, 
c 


EL, for each of the discrete horseshoe vortices. 


The valves of the increments of the lateral lift distributions A CE 


dy 


on the fuselage due to all of the horseshoe vortices and their images 
obtained from equation (5) are superimposed to get the complete lateral 


lift distribution dbp/dy. 
of horseshoe vortices and imsges 


The increment of total lift for each pair 


Ale 


and the longitudinal center of pressure of the total lift 


acts at the bound-vortex axis, 


Le is 


obtained simply by dividing the sum of the moments of the incremental 


total lifts ALr by the total lift 


Le (the sum of the increments ALp). 


The results of the calculation for the lateral distribution and 
longitudinal center of pressure carried out in the appendix are pre- 


sented in figures 10 and 11, respectively. 
contributed by each pair of vortices and images 


ment of total lift Alp 


Figure 11 shows the incre- 


and it can be seen that the contribution to the total lift Le due to 
the outboard part of the wing (the contributions farthest aft in the 


side view) are small compared to the 
of the wing. 


contribution of the inboard part 


The induced lift on the fuselage has also been calculated 


by using the wing-alone spanwise lift distribution shown in figure 9; 


ee 
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that is, the effect of the fuselage on the wing is neglected. It can 
be seen in figures 10 and 11 that, although there is about a lO-percent 
increase in the magnitude of dlp /dy and Lp due to the effect of the 


fuselage on the wing, the lateral and longitudinal load distribution on 
the fuselage is practically unaffected. The correction for finite length 
has not been included in the calculations. Its only effect would have 
been to decrease the magnitude of the lift, but it would not alter the 
lateral or longitudinal distribution. 


The total Lift and moment on the fuselage can be obtained by adding 
the components of lift and moment due to the induced lift to the com- 
ponents of lift and moment due to the local angle of attack of the fuse- 
lage (which can be calculated as shown in ref. 2). The lift and moment 
on the part of the wing outboard of the fuselage wall calculated by the 
method of reference 6 or the method of this paper (described in a later 
section) may then be added to the lift and moment on the fuselage to 
get the total lift and moment on the combination. 


Outline of Method for Calculating the Lift on the Wing in 
the Presence of the Fuselage 


The Lift distribution on the wing will be calculated by equating 
the downwash angle induced at the three-quarter-chord line by the horse- 
shoe vortices centered on the quarter-chord line (see fig. 5) to the 
local angle of attack on the three-quarter-chord line at several points 
along the span. Since the boundary conditions on the fuselage are not 
completely satisfied by the vortex-image system of figure 5 in the 
region near the bound vortex, it is necessary to resort to certain l 
approximations in calculating the downwash. Calculation of the exact 
values of these downwash functions would require the calculation of the 
additional potential which involves a great deal of time and effort, 
but if such calculations were made, the values could be used directly 
in the present method and the restrictions suggested below would be 
eliminated. 


The approximations described in the following sections improve as 
the longitudinal distance from the downwash points to the bound vortex 
increases so that this method is considered to be best-suited for con- 
figurations having fairly small ratios of diameter to root chord, say 

traight wings with -Dlameter <j and for swept wings with 
kor EE ngs Root chord 5 x " 

_Diameter s l, Calculations of the aerodynamic loading on a wing and 
Root chord 

tip-tank combination, using a method corresponding to the one described 
herein with the downwash points located about 1 tip-tank diameter behind 


7th 
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ーー に 5), have been found to yield results in 


Tip chord 2 
good general agreement with experimental results. 


the bound vortex ( 


Downwash in region near the wing bound vortex.- Although the 
boundary conditions are satisfied completely only at infinity, the 
trailing legs of a single pair of horseshoe vortices and images also 
satisfy the boundary conditions on the cylinder in the plane perpen- 
dicular to the cylinder axis which passes through the bound legs of the 
horseshoe vortices. In addition, the boundary conditions on the cylin- 
der are satisfied completely by the semi-infinite vortex-image system 
everywhere in the plane of the horseshoe vortex. Since the boundary 
conditions on the cylinder are satisfied exactly at the points noted, 
and are partly satisfied everywhere else, it will be assumed that the 
downwash due to the trailing vortices may be calculated approximately, 
at least in the plane of the horseshoe vortex, without introducing any 
correction factor. 


It will be necessary, however, to use a correction factor in calcu- 
lating the downwash due to the real and image bound vortices which have 
the greatest tendency to violate the boundary condition of zero velocity 
normal to the surface of the cylinder. It will be assumed that, along 
a line parallel to the bound vortex, the effect of the cylinder on the 
vertical flow induced by the bound vortex will be the same as its effect 
on a two-dimensional uniform rectilinear flow, so that the downwash 
induced by the bound vortices on that line is increased by the fac- 
tor 1+ =. 


2 
Ya 


The downwash angle at a point at y=ya and x= x, (shown in 
fig. 4) is then written | 


2 
a = T nya 2 Ya Fn( xy O) + Gn(x,y,0} i + a (18) 
a 


where the downwash factor due to the trailing vortices F, and the 
downwash factor due to the bound vortices Gp may be calculated by 
the Biot-Savart law. 
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The function F4 is the sum of four terms, each having the form 


(the term corresponding to the trailing legs of the horseshoe vortex 
centered at y = yn) and the function Gp is the sum of four terms, 


each having the form 


Yn - Ya + 1 yn - Ya - 1 


(Sn - 94 + 2)? + (x, - x)? (3a - Fa - 3)? + (ža - 3a) 


(the term corresponding to the bound leg of the horseshoe vortex 
centered at y = Yu). 


Scheme for calculating spanwise lift distribution on the wing.- The 
lift distribution on the wing is calculated by equating the downwash 
angle on the wing induced by the vortex-imsge system to the local angle 
of attack on the three-quarter-chord line of the wing at n points as 
in reference 6. Thus equation (18) is written for each of the n sta- 
tions, and the simultaneous equations can then be solved for yp at 
n points on the span. The + marks in figure 5 indicate the points 
where the downwash is equated to the local angle of attack. 


The effective angle of attack a, 1s equal to the geometric angle 
of attack on the wing og plus the angle of attack induced by the 


fuselage 


= Cg + ap’ ; l . (19) 
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and 


5 

| 

& 
UN 

+ 
Pol ® 

[ 

A 


a? (20) 


ty 


where ap is the geometric angle of attack of the fuselage and the 
factor 1 + s takes into account the increase in the vertical velocity 


of the free stream in the neighborhood of the fuselage as calculated by 
assuming the fuselage to be an infinite cylinder in a two-dimensional 
uniform rectilinear flow of magnitude Vap. 


The effect of finite fuselage length (referred to in ref. 6 as the 
"inflow effect") must be included separately by multiplying y at each 
spanwise station by the factor 1+ 25 where 5 is the ratio of the 
local increment of longitudinal velocity due to the fuselage to the free- 
stream velocity. As shown in reference 6, the factor 1+ 25 is used 
to account approximately for the small increase in dynamic pressure in 
the flow over the wing due to the increase in the local Longitudinal 
velocity near the surface of a fuselage of finite length. 


Outline of Method for Calculating the Downwash 


The downwash calculations made in this section are for the case 
when the wing is at an angle of incidence with the fuselage at zero 
geometric angle of attack. The accuracy of the calculated results must 
be verified by experimental results; however, by comparison with results 
of downwash calculations for the wing alone, the results of the calcu- 
lations for the downwash in back of the wing-fuselage combination may 
serve to give useful information regarding the effect of the wing-fuselage 
interference. 


As was noted previously, the boundary conditions on the wing-fuselage 
combination are completely satisfied by the vortex-image system infinitely 
far behind the wing. At a great distance behind the wing, therefore, the 
theoretical value of the downwash can be calculated accurately for a 
given spanwise lift distribution. This calculation is simple, and the 
results may be useful for certain applications. However, the approxi- 
mations made in the preceding section for calculating the downwash at 
the wing three-quarter-chord line may give more accurate results for 
the downwash calculated in the plane of the wing in the region nearer 
the bound vortex, 


——— 
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If the lift distribution on the wing is known, say from the method 
of reference 6, or the method outlined in this paper, the downwash angle 
in the plane x = o may be calculated by adding the downwash due to 
each pair of trailing vortices and their images. Thus, from the Biot- 
Savart Law and figure 4, the downwash angle a, at the point y = Ya» 


Z = Za, and x= co is given as 


1 
Ta Sar 2 nn (i 


where 


Xn ~ Yg + l Yn - Ye - 1 


x pa 2 = m = 2 - 2 
(Fn ey bb 1) + っ. ぞ due 9 1) + Za 


Yn E Ya +1 Yn + ya - 1 


— mo 
- = ニー + L ¥, - x =- 1l 
i n a Jn m 
z 2 2  _ - = | ご 
Ja ご ミー + 1 + A (Fa - se ー + z ぞ 
Yn X Jn 
= の =D 
Şa +5 +1 J, E- 
Yn Yn 
(22) 
- a2 = 2 の 2 2 
Ya エニー エエ ユエ) +5 Ya Pag tag l] + Ze 
Jn Jn 


and 7, = (ce;/8), is the local lift at station typ. 
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The representation of the wake by discrete vortices is satisfactory 
if the downwash is calculated at points halfway between the trailing 
Legs (Ya = Yn). The values calculated at these points may be faired 


to obtain a continuous spanwise distribution of downwash angle. 
CONCIUDING REMARKS 


For the purpose of calculating the longitudinal loading on the fuse- 
lage in subsonic flow, a midwing wing-fuselage combination has been repre- 
sented by a system of discrete horseshoe vortices and images. Using this 
simplified mathematical model, a method for calculating the lift on the 
fuselage induced by the wing is presented and illustrated by a numerical 
example. This "induced 1lift™ can be added to the lift on the part of 
the wing outboard of the fuselage and the lift on the fuselage due to 
the upwash induced by the wing to get the total loading on the wing- 
fuselage combination with an axisymmetrical fuselage and a wing of 
arbitrary plan form. 


ln addition to the method for calculating the induced lift, which 
is theoretically rigorous, methods for calculating the downwash far 
behind the wing and for calculating the spanwise lift distribution on 
the wing in the presence of the fuselage are outlined. In calculating 
the spanwise lift distribution on the wing, approximations are made to 
account for the effect of the "additional potential" so that the method 
is not rigorous. However, the effect of the additional potential may 
easily be incorporated into the method when it is calculated. 


Langley Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va. 
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APPENDIX 


ILLUSTRATIVE EXAMPLE 


A numerical example is given to illustrate the method presented in 
the body of the paper for calcuiating the induced lift on the fuselage. 


Geometric characteristics of the midwing configuration used in the 


illustrative example. - The 


is essentially the same as 


Aspect ratio . 
Taper ratio . 
Sweepback, deg 
a* 
g* 


plan view of the wing-fuselage combination 
shown in figure 5. The geometric data are 


e 8 e 8 L] 
[| [| a 8 s 
. L] 6 8 s 
. 8 e $ . 
a a $ Ld a 
LI a 8 ® ? 
6 a | | 8 LI 
|] a 8 8 合 
6 . 8 LI . 
a a e ¢ | 
e. . a a 8 
e. * » . LI 
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Spanwise loading on the wing in the presence of the fuselage.- The 
spanwise Lift distribution on the wing in the presence of the Puselage 
tabulated below and shown in figure 9 is the same as the one obtained 
in reference 6 for a high midwing configuration. It is, however, quali- 
tatively correct for the pure midwing case and is used to illustrate 
the procedure for obtaining the induced lift on the fuselage. The lift 
distribution on the wing in the presence of the fuselage plotted in fig- 
ure 9 is tabulated as follows: 


1 
2 
3 
ip 
5 
6 
7 
6 
9 


Longitudinal distribution of the induced lift.- Equation (6) which 
gives the lift on the fuselage due to a single pair of horseshoe vortices 
at (Xn' tyn’) and their images can be written 


Tupa € 


22 ORAS NACA RM L52Je Ta 


SR MEGAN. a 3 (A1) 
n y. 2 gta 


sey I の 


n Yu? a gt 


The longitudinal center of pressure referred to the intersection of the 
quarter-chord line and the fuselage axis of symmetry (gee fig. 5) is 


2 


n Yu!” € g1 


3 (=) 2 
C n 1 2 2 


C 
n y TE 


Xop (43) 


The computed values of Lp/qS and Xop' are obtained from the values 
presented in the following table: 


1 
2 
3 
h. 
2 
6 
T 
8 
9 
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n 
O 
+ 
Q 
UY 
no 


: n Jn = a! 
Zep = 
(9 + 
l C /n 2 
X: Bt 
= 2.51 
ce i 
The values of s*(— i are plotted as vectors in fig- 
Č /n y, z gi? 


ure 11. The location of Xop' is also shown. 


Lateral distribution of the induced lift.- Equation (5) which gives 
the lateral distribution of the induced lift on the fuselage due to a 
single pair of horseshoe vortices at (=, ' »iYn') and their images can be 
written 


(xa' - s)? - 1 (yu? +s") - 1 


(Ah) 


^i) Do. (で) に E 2(yn' - s") - yy" - 8 りや - y= eet 2(yn' + s') NI -y2 


The total lift at a station ty' on the fuselage is then 


a -25 En n71 im - -2)h - p* iuro + mW - y!? 
qe NE で NS on = i 2. (Ya' A s')? -1 


(A5) 
! 


dl. 


The value of ーー at ty! - 0,25 is calculated from the values in the following table (the 
q 


He 


lateral lift at any station ty! can be calculated in a similar manner): 


CC7 


NO 00-93 Oo FW NO FJ 


at y' = 0,25. The complete lateral lift distribution is shown in figure 10. 
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Figure 1.- Combinstion of horseshoe vortex and infinite cylinder. 
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Figure 2.- Illustration showing superposition of infinite trailing 
vortices on cylinder and horseshoe vortex system. 


ef 


28 


AA) NACA RM L52J2'T& 


JV 
ンー へ ヽ 

も に 

5 
ict 


| | 
| | 


D wl | | 
| 
y 


X 


Figure 3.- Superposition of two unequal pairs of symmetrical horseshoe 
vortices to form & discrete horseshoe vortex-image system. 
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Figure 4.- Cylinder with a pair of symmetrically disposed horseshoe 
vortices and their images. 
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Sweepback angle, 45°; 


a* = 0.10. 


taper ratio, 0.45; A = 8; 


Figure 5.- Cylinder-wing combination. 
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(b) Infinite-vortex—ellipsoid configuration. 


Figure 6.- Illustrations of sphere with infinite vortex end ellipsoid 
with infinite vortex. 
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Figure 7.- Effect of slenderness ratio on the maximum velocity of flow 


past an ellipsoid of revolution, 
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Figure 8.- Effect of slenderness ratio on the lift carried by an ellipsoid 
of revolution intersected by an infinite vortex having the strength Ih. 
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Figure 9。- Lateral lift distribution on wing alone and on a wing-fuselage 
combination. Sweepback angle, 159; taper ratio , 0.45: A = 8; at = 0.10. 
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Figure 10.- Lateral distribution of induced lift on the fuselage in the 
presence of the wing. Sweepback angle, 45°; taper ratio, 0.45; A = 8; 
&* = 0,10. . 
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Figure 11.- Longitudinal distribution of induced lift on fuselage in the 
presence of the wing. Sweepback angle, 45°; taper ratio, 0.45; A = 8; 
a* = 0,10. 
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